Pion correlations in hydro-inspired models with resonances 

WojciECH Florkowski, Wojciech Broniowski 

Institute of Physics, Swigtokrzyska Academy, 
ul. Swigtokrzyska 15, 25-406 Kielce, Poland, 
and 

Institute of Nuclear Physics, Polish Academy of Sciences, 
ul. Radzikowskiego 152, 31-342 Krakow, Poland 

AND 

Adam Kisiel, Jan Pluta 
Faculty of Physics, Warsaw University of Technology, 00-661 Warsaw, Poland 



The effects of shape of the freeze-out hypersurface and resonance de- 
cays on the pion correlation functions in relativistic heavy-ion collisions are 
studied with help of the hydro-inspired models with single freeze-out. The 
heavy-ion Monte-Carlo generator THERMINATOR is used to generate hadronic 
events describing production of particles from a thermalized and expanding 
source. We find that the short-lived resonances increase the pionic HBT 
radii by about 1 fm. We also find that the pion HBT data from RHIC 
are fully compatible with the single freeze-out scenario provided a special 
choice of the freeze-out hypersurface is made. 

PACS numbers: 25.75.-q, 25.75.Dw, 25.75.Ld 

1. Introduction 

In this lecture we discuss the pion correlation functions ^ 121 131 13 El El 

Ej ) obtained from the hydro- inspired statistical model of hadronization 
implemented in THERMINATOR 9 . Our theoretical predictions for relativistic 
Au -|- Au collisions at RHIC top energies are compared with the STAR data 
jlOj and a good agreement is found for a freeze-out hypersurface with the 
transverse radius decreasing with time. Our presentation is based on the 
recent paper on femtoscopy in hydro-inspired models with resonances |11| . 

* Presented by WF at XLVI Cracow School of Theoretical Physics, 27.05.06 - 5.06.06, 
Zakopane. 
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The hydro-inspired models use concepts borrowed from relativistic hy- 
drodynamics but they do not include the complete time evolution of the 
system. Such models help us to verify the idea that matter, just before 
the kinetic freeze-out is locally thermalized and exhibits collective behavior 
[121 113j . The observables are expressed in terms of thermal (Bose-Einstein, 
Fermi-Dirac) distributions supplemented with the collective expansion of 
the system and decays of resonances. 

Following Refs. |14| ll5j. in our approach we assume one universal freeze- 
out for all processes; inelastic and elastic processes cease at the same time, 
also emission of strange and ordinary hadrons occurs at the same moment. 
This is a simplifying but very fruitful assumption. ^ It has been shown in a 
series of papers |[17i tlSj llSl t20|l |2lj that this approach gives good description 
of particle yields, transverse-momentum spectra, pion invariant-mass dis- 
tributions, balance functions, azimuthal asymmetry, and transverse energy. 
The single freeze-out model is also consistent with the sudden hadronization 
(explosion) scenario at RHIC proposed in Ref. 22.. 

In this lecture the results obtained with the Monte-Carlo implementation 
of the single-freeze-out model are presented 9 . This method allows us to 
use a two-particle method to extract the correlation functions. Moreover, 
the Monte-Carlo approach forms a very convenient platform for the inclusion 
of Coulomb effects. 



For boost-invariant and cylindrically symmetric models the freeze-out 
hypersurface is defined by the freeze-out curve in Minkowski space t — p. 
Such a curve is obtained by the projection of the freeze-out hypersurface on 
the plane = Q |121 113j . In our analysis we consider the Cracow model and 
the (generalized) blast wave-model with resonance decays. The freeze-out 
hypersurface and the flow for the Cracow model have the following form 



2. Freeze-out hypersurface and flow 



.2 2 

t -r. 



+ , T = const. 



(1) 





^ A recent study of Nonaka and Bass, Fig. 16 in Ref. |16| . shows that rescattering 
effects have httle impact on the particle spectra. Rescattering effects on particle 
correlations could be studied with a suitable "afterburner" for THERMINATOR. 
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Fig. 1. Various parameterizations of the freeze-out hypersurface. The curves show 
the dependence of time t on the radial distance p - 
four models considered. 
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The blast wave model uses the following parameterization 

(r + ap)^, r, a = const, 
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v± = const. 



In Fig. ^we show the freeze-out curves for the four considered models; the 
Cracow model and the blast-wave model with three different options for the 
parameter a (a = 0.5, a = 0, a = —0.5). All these forms describe properly 
the transverse-momentum spectra, however, as we shall see below, they lead 
to different predictions for the pion correlation functions. 

It is important to emphasize that the Cracow and blast-wave models 
are treated by us on the same footing, i.e., the effects of the resonance 
decays are included in the all considered models (our calculations take into 
account all well established resonances, 381 particle types with the total of 
1872 decay modes, we note that this input is the same for both THERMINATOR 
and SHARE 23 ). The only important difference between the models resides 
in the definition of the freeze-out hypersurface, see Eqs. © - ©• When the 
freeze-out hypersurface and the flow are defined, the Cooper-Frye formula 
is used to construct the emission function S{x,p), which is later used to 
calculate various physical observables. Details of the construction of the 
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emission function including the effects of the resonance decays are given in 
Ref. HI]. 

3. Correlation functions 

3.1. Basic definitions 
The measured correlation function is given by the formula 

g Pi,P2 = ..r (3) 

Wi{pi)Wi{p2) 

where Wi and W2 are one- and two-particle distributions. The model cal- 
culations relate the correlation function to the emission function, 

^) ^ ! d''xiS{xi,pi)d^X2S{x2,P2)\^{k*,n? 

j d'^XiS{xi,pi) j d'^X2S{x2,P2) 

where \^{k* ,r*)\'^ is the squared wave function of the pion pair, q is the 
relative momentum of the two particles forming a pair, k is the average 
momentum, k* is the momentum of the first particle in the pair rest frame, 
and r* is the relative separation in this frame. 

3.2. Monte-Carlo Method 

In our calculation, which is based on the Monte-Carlo method, the in- 
tegration in Eq. Q) is replaced by the summation over particles and pairs 
of particles 

E E 5^{q-Pi+m^{k-\{p, +pj))\m{k*,r*)? 
Ciq, k) = ^ , (5) 

EE'JA(g - pi +Pj)<5A(fe - ^{pi +Pj)) 

i j 

where d^ip) is a box-like function defined by the expression 

1 if \Px\ < f , \Py\ < \Pz\ < t 

Sa{p) = { (6) 
otherwise. 

In the numerical calculations we use the bin value A = 5 MeV. The ranges 
of k are taken from STAR experiment: 0.15 - 0.25 GeV, 0.25 - 0.35 GeV, 
0.35 - 0.45 GeV, and 0.45 - 0.60 GeV. 

For each pair, considered in Eq. the following transformations are 
made: Firstly, the pair is boosted from the laboratory frame to the longi- 
tudinal co-moving system (LCMS) using the Bertsch-Pratt decomposition 
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into out-, side-, and long- components. Secondly, the pair is boosted from 
LCMS to the pair rest frame (PRF). In this way the correlation function 
becomes a histogram of the squares of the wave function calculated for each 
pair in PRF where it is defined, but tabulated in LCMS. 



3.3. Wave functions 

We consider two options for the wave function: i) the simplest wave 
function includes symmetrization over the two identical pions but neglects 
any dynamical interactions. 



1 



Jk*r* 



q,Q = (e- ' +e" 



■ik*r* 



(7) 



and ii) the Coulomb interaction is included, and the wave function has the 
form 



-ik*r* 



F(-Z77,l,ie+) 



(8) 



where 5c is the Coulomb phase shift, Ac is the Coulomb penetration factor 
(sometimes called the Gamow factor), = k*r* it k*r* = k*r*{\ it cos0*), 
rj = {k*a)~^ with a being the Bohr radius of the pair, and F is the confluent 
hypergeometric function. The angle between k* and r* is denoted by 6* . 

If the simple wave function ((T)) is used, the 3D correlation function is 
fitted with the standard gaussian formula 



C = 1 Aexp -i?out(^-L)9out - ^side(^-L)9side " Aong(^-L)9j 



ong 



(9) 



On the other hand, when the Coulomb wave function (jSJ is used, the 3D 
correlation function is fitted with the Bowler-Sinyukov formula [241 125j 



C{q, k) = {l-\)+ AKcoul(ginv) 1 + exp (--Rout^out 



(10) 



where i^coui('?inv) is the squared Coulomb wave function integrated over a 
static gaussian source. We use, following the STAR procedure ^0], the static 
gaussian source characterized by the widths of 5 fm in all three directions 
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Fig. 2. Transverse-mass spectra at mid-rapidity of pions (open circles) and kaons 
(open squares) for the Blast- Wave model with resonances and a = —0.5, T = 165.6 
MeV, HB = 28.5 MeV, r = 8.55 fm, p^ax = 8.92 fm, and v± = 0.311 c. The data 
points (stars for pions, triangles for kaons) come from the STAR collaboration |10|. 



In our calculations the parameters of each model are fixed by fitting the 
single-particle p^-spectra of pions and kaons to the experimental data. An 
example of such a fit is shown in Fig. [2 The values of the parameters are 
given in the captions of Figs. |21 and 0] where our main results concerning 
the intercept A and the HBT radii -Routi ^sidej and -Riong are presented 
as functions of the transverse momentum of the pion pair. The squares 
correspond to the full calculation with resonances, the down-triangles show 
the results obtained in the calculation without resonances, the up-triangles 
show the results obtained with resonances and with the Coulomb-aware 
fit made according to the Bowler-Sinyukov formalism [241 I25j . while the 
circles show the data of the STAR collaboration from Ref. ^U]. One can 
observe that the inclusion of resonances increases the radii by about 1 fm. 
This is expected, since the resonances travel some distance from their place 
of birth on the freeze-out hypersurface before they decay into pions. The 
typical scale is set by the resonance life-time which is about 1 fm. We 
observe a decrease of the radii with kx which is a known qualitative effect 
of the presence of the flow correlating momenta with emission points. We 
also note that the effect of resonance decays is larger at small kx- 

The model values of the intercept A shown in Figs. I3I4I are too large 
compared to the data, which simply reflects the fact that we do not take 



4. Results 
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Fig. 3. (Left) The results for Cracow model: A and the HBT radii i?out, -Rsidc, and 
i?iong shown as functions of the transverse momentum of the pion pair. The squares 
show the full calculation with resonances, down-triangles is the same without res- 
onances, the up-triangles show the calculation with resonances and the Coulomb 
corrections made according to the Bowler-Sinyukov method, while the circles show 
the data of the STAR collaboration for ^Jsnn = 200 GeV ^Uj- The lines are 
drawn to guide the eye. We note that the inclusion of resonances increases the 
radii by about 1 fm. The model parameters are: T — 165.6 MeV, /is = 28.5 MeV, 
T = 10.55 fm, and pmax = 7.53 fm. (Right) Same as the left panel but for the 
Blast-Wave model with resonances and a — 0.5. The model parameters in this 
case are: T = 165.6 MeV, ^lB = 28.5 MeV, r = 9.91 fm, p^ax = 7.43 fm, and 
vi_ 0.407 
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Fig. 4. (Left) Results for the Blast- Wave model with resonances and a — 0. The 
model parameters are: T — 165.6 MeV, = 28.5 MeV, r = 8.17 fm, pmax = 
8.21 fm, and v± — 0.341. (Right) Results for the Blast- Wave model with resonances 
and a ~ —0.5. The model parameters are: T — 165.6 MeV, /is — 28.5 MeV, 
T = 8.55 fm, Pmax = 8.92 fm, and v± = 0.311 . This is the model that produces 
best agreement out of the four models considered. With the same values of the 
parameters the model reproduces the transverse-mass spectra. Open squares in 
the top panel show results for the analysis of A including all weak-decay products. 

into account the effect of secondary pions coming from the weak decays, as 
well as the contamination of the pion sample by misidentified particles in 
the experiment. The inclusion of the weakly-decaying particles is presented 
in Fig. |1] (open squares). We see a dramatic drop of the A parameter, as 
expected. The remaining discrepancy between model values and experiment 
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can be attributed to the misidentification of particles in the experiment. 

We emphasize that the results shown in Figs. 01 and 0] are obtained by 
fitting the 3-dimensional two-particle correlation function. This procedure 
reveals non-gaussian features of the correlation functions. In particular, the 
pairs where one of the pions comes from the oo decay produce long-range 
tails caused by the long lifetime of the uj meson. Such features are discussed 
in more detail in Ref. 

5. Conclusions 

Our results show that simultaneous description of the transverse - mo- 
mentum spectra and the correlation radii is possible in the hydro-inspired 
models if special choice of the freeze-out hypersurface is made. We have 
found that the data favor the freeze-out geometry where the transverse size 
decreases with time. 

Our approach is as close as possible to the experimental treatment of 
the correlations; we use the two-particle method and include the Coulomb 
corrections. The role of the resonances has been analyzed in detail. Some 
earlier expectations concerning lowering of the intercept and the role of 
the omega meson have been confirmed. On the other hand, in contrast to 
earlier studies, we have found that the strong decays of resonances increase 
the radii by about 1 fm. 
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